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1 Introduction 

We introduce a new problem, CSPLib problem number 50, to generate all degree 
sequences that have a corresponding diamond-free graph with secondary prop- 
erties. This problem arises naturally from a problem in mathematics to do with 
balanced incomplete block designs; we devote a section of this paper to this. The 
problem itself is challenging with respect to computational effort arising from 
the large number of symmetries within the models. We introduce two models 
for this problem. The second model is an improvement on the first, and this 
improvement largely consists of breaking the problem into two stages, the first 
stage producing graphical degree sequences that satisfy arithmetic constraints 
and the second part testing that there exists a graph with that degree sequence 
that is diamond-free. We now present the problem in detail and then give moti- 
vation for it. Two models are then presented, along with a listing of solutions. 
We then conclude and suggest further work that might be done. 

2 Problem Definition 

Given a simple undirected graph G = (V,i£), V is the set of vertices and E 
the set of undirected edges. The edge {u, v} G E if and only if vertex u is 
adjacent to vertex v in G. The graph is simple in that there are no loop edges, 
i.e. \/ v ^v [{ v , v } £ E\. Each vertex v in V has a degree 5(v) = |{{v, w} : {v, w} G 
E}\, i.e. the number of edges incident on that vertex. A diamond is a set of four 
vertices in V such that there are at least five edges between those vertices (see 
Figure [l] for an example of a diamond). Conversely, a graph is diamond-free if it 
has no diamond as a subgraph, i.e. for every set of four vertices the number of 
edges between those vertices is at most four. 

In our problem we have additional properties required of the degree sequences 
of the graphs, in particular that the degree of each vertex is greater than zero 
(i.e. isolated vertices are disallowed), the degree of each vertex is divisible by 3, 
and the sum of the degrees is divisible by 12 (i.e. \E\ is divisible by 6). 

The problem is then for a given value of n, such that \V\ = n, produce 
all unique degree sequences 6(1) > 5(2) > ... > S(n) such that there exists 
a diamond- free graph with that degree sequence, each degree is non-zero and 
divisible by 3, and the number of edges is divisible by 6. 
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Fig. 1. A simple diamond graph of four vertices and five edges. 

In Figure [2] we give the unique degree sequence for n = 8 and an adjacency 
matrix and simple graph that both corresponds to that sequence and represents 
a diamond-free graph. 



33333333 

111 
00000 111 
000 1100 1 
10 10 10 
1 10 10 
110 10 
110 10000 
11100000 




Fig. 2. A degree sequence for n = 8 with the corresponding adjacency matrix 
and graph that is diamond- free. 



3 Motivation 

The problem is a byproduct of attempting to classify partial linear spaces that 
can be produced during the execution of an extension of Stinson's hillwalking 
algorithm for block designs with block size 4. First we need some definitions (see 

HD- 

Definition 1. A Balanced Incomplete Block Design (BIBB) is a pair (V,B) 
where V is a set of n points and B a collection of subsets of V (blocks) such 
that each element of V is contained in exactly r blocks and every 2- subset of V 
is contained in exactly A blocks. 

Variations on BIBDs include Pairwise Balanced Designs (PBDs) in which blocks 
can have different sizes, and linear spaces which are PBDs in which every block 
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Algorithm 1 Algorithm to generate an STS on n points 
n <(— number of points 
LivePairs <(— {(i, j) : < i < j < n} 
Blocks <(— empty set 
while LivePairs not empty do 

choose (x, y) and (y, z) from LivePairs 

remove (x, y) and (y, z) from LivePairs 

add (x, y, z) to Blocks 

if (y, z) is in LivePairs then 
remove (y, z) from LivePairs 

else 

remove existing block containing (y, z), (w, y, z) 
add (w, z) to LivePairs 
end if 
end while 



has size at least 2. It is usual to refer to the blocks of a linear space as a line. 
A partial linear space is a set of lines in which every pair appears in at most A 
blocks. Here we refer to a BIBD with A = 1 as a block design and to a partial 
linear space with A = 1, having lines of size z, where i > 3 and s^ > as a 
3 S3 4 S4 . . . structure. For example, a block design on 7 points with block size 3 is 
given by the following set of blocks: 

(0, 1, 2), (0, 5, 6), (0, 3, 4), (1, 4, 5), (2, 3, 5), (2, 6, 4), (1, 3, 6) 

and a 3 4 4 : structure by the set of blocks 

(0,1, 2, 3), (0,4, 5), (0,6, 7), (1,4, 6), (1,5, 7) 

Note that in the latter case we do not list the lines of size 2. Block designs with 
block size 3 are known as Steiner Triple Systems (STSs). These exist for all n 
for which n = 1, 3 (mod 6) [4]. For example the block design given above is the 
unique STS of order 7 (STS (7)). Similarly block designs with block size 4 always 
exist whenever n = 1,4 (mod 12). 

Algorithm [I] allows us to generate an STS for any n and is due to Stinson 
[5 . This algorithm always works, i.e. it never fails to terminate due to reaching 
a point where the STS is not created and there are no suitable pairs (x, y) and 
(x,z). 

A natural extension to this algorithm, for the case where block size is 4, is 
proposed in Algorithm [2] This algorithm does not always work. It is possible for 
the algorithm to fail to terminate due to reaching a point where the block design 
is not created and there are no suitable overlapping triples (x, y, z) and (x, y, w) 
in Live Triples. For this reason, we replace the condition on the while loop by 

while LiveTriples not empty and overlapping triples exist. 

Now the algorithm terminates, but rather than always producing a block design, 
either produces a block design, or a 4 S4 structure, for which the complement has 
no overlapping triples. I.e. the complement graph is diamond-free. 
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Algorithm 2 Algorithm to generate a block design with block size 4 on n points 
n <(— number of points 

LiveTriples <(— {(i, j, k) : < i < j < k < n} 

Blocks <(— empty set 

while LiveTriples not empty do 

choose (x,y,z) and (x,y,w) from LiveTriples 

remove (x,y,z) and (x,y,w) from LivePairs 

add (x,y,z,w) to Blocks 

if (y, z, w) is in LiveTriples then 
remove (y, z, w) from LiveTriples 

else 

remove existing block containing (y, z, w), (u, y, z, w) 
add (u,y,z) and (u, z,w) to LiveTriples 
end if 
end while 



When n = 13 the algorithm either produces a block design, or a structure 
whose complement graph consists of 4 non-intersecting triangles. 

The next open problem therefore is for n = 16. If the algorithm does not 
produce a block design, what is the nature of the structure it does produce? 
To do this, we need to classify the 4 r4 structures whose complement graph is 
diamond- free. 

The cases for which the 4 S4 structure has at least 2 points that are in the 
maximum number of blocks (5) are fairly straightforward. (There are fewer cases 
as this number increases.). However if the number of such points is or 1, 
there is a large number of sub-cases to consider. The problem is simplified if 
we can dismiss potential 4 S4 structures because the degree sequences of their 
complements can not be associated with a diamond-free graph. This leads us to 
the problem outlined in this report: to classify the degree sequences of diamond- 
free graphs of order 15 and 16. Note that each point that is not in 5 blocks is 
either in no blocks or is in blocks with in some number of points, where that 
number is divisible by 3. Thus for every point there is a vertex in the complement 
graph whose degree is also divisible by 3. In addition, since the number of pairs 
in both a block design on 16 points and a 4 S4 structure are divisible by 6, the 
number of edges in the complement graph must be divisible by 6. 

4 Constraint Models for Diamond-free Degree Sequences 

We present two constraint models for the diamond-free degree sequence problem. 
The first model we call model A, the second model B. In many respects the 
two models are very similar but what is different is how we solve them. In the 
subsequent descriptions we assume that we have as input the integer n, where 
\V\ = n and vertex i E V. All the constraint models were implemented using the 
choco toolkit [2]. Further we assume that a variable x has a domain of values 
dom(x). 
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4.1 Model A 

Model A is based on the adjacency matrix model of a graph. We have a 0/1 
constrained integer variable for each edge in the graph such that = 1 
{i, j} G E. In addition we have constrained integer variables degi to deg n corre- 
sponding to the degrees of each vertex, such that 

V iG [i.. n ] dom(degi) = [3 .. n - 1] (1) 
We then have constraints to ensure that the graph is simple: 

V;e[i..n]Vj G [;.. n ] [Aij = A^] (2) 

Vi G[ l..n][Ai = 0] (3) 

Constraints are then required to ensure that the graph is diamond-free: 

V{i,j,fc,z}Gv[^i,j + A i,k + A»,z + A jik + A i; , + A ki i < 4] (4) 
Finally we have constraints on the degree sequence: 

j=n 

V iG [i .. n][de& = (5) 

Vi G [i .. n -i][de^i > de^+i] (6) 
V iG [i .. n ][de^ mod 3 = 0] (7) 

a = ^ de#i (8) 
a mod 12 = (9) 
The vertex degree variables degi to deg n are the decision variables. 



4.2 Model B 

Model B is essentially model A broken into two parts, each part solved separately. 
The first part of the problem is to produce a graphical degree sequence that 
meets the arithmetic constraints. The second part is to determine if there exists 
a diamond-free graph with that degree sequence. Therefore solving proceeds as 
follows. 

1. Generate the next degree sequence tt = di, c^, d n that meets the arith- 
metical constraints. If no more degree sequences exist then terminate the 
process. 

2. If the degree sequence tt is not graphical return to step 1. 

3. Determine if there is a diamond- free graph with the degree sequence tt. 

4. Return to step 1. 



6 Diamond-free Degree Sequences 



The first part of model B is is then as follows. Integer variables degi to deg n 
correspond to the degrees of each vertex and we satisfy constraints (1), (6), (7), 
(8) and (9) to generate a degree sequence. 

Each valid degree sequence produced is then tested to determine if it is graph- 
ical using the Havel-Hakimi algorithm [3] . If the degree sequence is graphical we 
create an adjacency matrix as in (2) and (3) and post the constraints (4) and 
(5) (diamond free with given degree sequence) where the variables degi to deg n 
have been instantiated. Finally we are in a position to post static symmetry 
breaking constraints. If we are producing a graph and degi = degj then these 
two vertices are interchangeable. Consequently we can insist that row i in the 
adjacency matrix is lexicographically less than or equal to row j. Therefore we 
post the following constraints: 

V*e[i .. n -i][degi = deg i+1 A { < A i+1 ] (10) 

where ■< means lexicographically less than or equal. In this second stage of 
solving the variables to A njU are the decision variables. 

5 Solutions 

Our results are tabulated in Table [I] (at end of report) for 8 < n < 16. All our 
results are produced using model B run on a machine with 8 Intel Zeon E5420 
processors running at 2.50 GHz, 32Gb of RAM, with version 5.2 of linux. The 
longest run time was for n = 16 taking about 5 minutes cpu time. Included in 
Table [I] is the cpu time in seconds to generate all degree sequences for a given 
value of n. 

All our results were verified. For each degree sequence the corresponding adja- 
cency matrix was saved to file and verified to correspond to a simple diamond- free 
graph that matched the degree sequence and satisfied the arithmetic constraints. 
The verification software did not use any of the constraint programming code. 

6 Conclusion 

We have presented a new problem, the generation of all degree sequences for 
diamond free graphs subject to arithmetic constraints. Two models have been 
presented, A and B. Model A is impractical, whereas model B is two stage and 
allows static symmetry breaking. 

There are two possible improvements. The first is to model A. We might add 
the lexicographical constraints, as used in model B, conditionally during search. 
The second improvement worthy of investigation is to employ a mixed integer 
programming solver for the second stage of model B. 

We are currently using the lists of feasible degree sequences for diamond-free 
graphs with 15 or 16 vertices to simplify our proofs for the classification of 4 S4 
structures with diamond-free complements, when the number of points in the 
maximum number of blocks is 1 or respectively. The degree sequence results 
for a smaller number of points will also help to simplify our existing proofs for 
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cases where more points are in the maximum number of blocks. Ultimately we 
would like to use our classification to modify the extension of Stinson's algorithm 
for block size 4 to ensure that a block design is always produced. 

In the more distant future, we would like to analyse the structures produced 
using our algorithm when n > 16. The next case is n = 25 and the corresponding 
diamond-free graphs would have up to 25 vertices. 
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n 


time 


deg 


ree sequence 




8 0.1 


33333333 


9 


0.1 


666333333 


10 


0.5 


6633333333 


11 


0.8 


63333333333 


12 


1.4 


333333333333 
666633333333 
666666666666 
966333333333 




13 


3.7 


6663333333333 
6666666333333 
6666666666633 
9633333333333 




14 


14.0 


66333333333333 
66666633333333 
66666666663333 
66666666666666 
93333333333333 
96666333333333 
99663333333333 
99933333333333 


15 107.7 


6 3 


3 


3333333333 


3 3 






6 6 


6 


6633333333 


3 3 






6 6 


6 


6666663333 


3 3 






6 6 


6 


6666666666 


3 3 






9 6 


6 


6333333333 


3 3 






9 6 


6 


6666633333 


3 3 






9 6 


6 


6666666663 


3 3 






9 9 


6 


3333333333 


3 3 






9 9 


6 


6666333333 


3 3 






9 9 


6 


6666666633 


3 3 






9 9 


9 


6663333333 


3 3 






9 9 


9 


9996666666 


6 6 






12 66333333333333 






12 12 12 333333333333 


16 339.8 


3 3 


3 


3333333333 


3 3 3 






6 6 


6 


6333333333 


3 3 3 






6 6 


6 


6666633333 


3 3 3 






6 6 


6 


6666666663 


3 3 3 






6 6 


6 


6666666666 


6 6 6 






9 6 


6 


3333333333 


3 3 3 






9 6 


6 


6666333333 


3 3 3 






9 6 


6 


6666666633 


3 3 3 






9 9 


3 


3333333333 


3 3 3 






9 9 


6 


6663333333 


3 3 3 






9 9 


6 


6666666333 


3 3 3 






9 9 


6 


6666666666 


6 3 3 






9 9 


9 


6633333333 


3 3 3 






9 9 


9 


6666666666 


3 3 3 






9 9 


9 


9333333333 


3 3 3 






9 9 


9 


9666666663 


3 3 3 






9 9 


9 


9666666666 


6 6 6 






9 9 


9 


9966666666 


6 6 3 






9 9 


9 


9996666666 


6 3 3 






12 633333333333333 






12 996333333333333 






12 12 66333333333333 






12 12 


93333333333333 



Table 1. Degree sequences, of length n, that meet the arithmetic constraints 
and have a simple diamond- free graph. Tabulated is n, cpu time in seconds to 
generate all sequences of length n and those sequences. 



